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Abstract 



We propose a simple approach for the radiative evolution of generic orbits around a Kerr black hole. For 
a scalar-field, we recover the standard results for the evolution of the energy E and the azimuthal angular 
momentum L z . In addition, our method provides a closed expression for the evolution of the Carter constant 
Q. 



Recently, the prospects of a direct observation of gravitational radiation led to a new interest in the 
old-standing problem of radiation back-reaction. The most promising radiation sources for the proposed 
gravitational- waves detector LIGO (as well as LISA [1]) are coalescing binary systems. [2] In such a binary 
system, energy and angular momentum are carried away by the emitted gravitational radiation. As a 
consequence, the system becomes more and more bound, resulting in a steady grow in the orbital frequency 
and energy-loss rate. Eventually, the binaries will coalless, releasing a strong burst of gravitational radiation. 
[2] The prospects of observing such bursts raises the need for a theoretical prediction of the evolution of such 
binary systems. If one of the binary objects is a black hole, then curvature effects could be important. This 
motivates one to understand the radiation-reaction in curved spacetime. 

Consider the situation of a compact object with mass /i orbiting a super-massive Kerr black hole with 
mass M » fi and angular momentum aM (we are using general-relativistic units C=G=1). This situation 
is particularly important for the proposed space-based gravitational- waves detector LISA, [1,3] but is also 
interesting in its own right. The problem of radiative evolution of orbits around a Kerr black hole was treated 
by several authors [4, 5], but so far with partial success only. The orbits in Kerr geometry are characterized 
by three constants of the motion: the energy E , the azimuthal angular momentum L z , and the Carter 
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constant Q . [6] Within the adiabatic approximation [3] that we use here, the problem of radiative evolution 
is reduced to that of calculating E , L z and Q , where the overdot denotes the evolution rate (in terms 
of external time t). So far, expressions have been given for E and L z , based on the rate that energy and 
angular-momentum are carried away to infinity and into the black hole [4,5,7]. However, Q is still unknown. 

[8] 

Although we arc primarily interested in the backreaction of gravitational radiation, it would be conve- 
nient to generalize the discussion and to consider electromagnetic and scalar-field radiation as well. We shall 
thus assume that the particle carries a "generalized charge" q , which may be either a scalar, electromagnetic, 
or gravitational charge (we demand q « M; in the gravitational case, q = [i ). In the discussion below 
we shall mainly refer to scalar and electromagnetic radiation (the explicit calculations in this Letter will be 
restricted to the scalar-field case) ; However, our approach can be extended to gravitational radiation as well 

Our main goal here is to propose a scheme for the calculation of Q (as well as E and L z ) for generic orbits 
in Kerr spacetime, by a direct calculation of the radiation-backreaction force - based on the pure retarded 
potential. An analogous direct scheme for radiative evolution of E and L z has already been proposed by 
Gal'tsov [4] , but that scheme is based on the so-called "radiative potential" (the "half- retarded- minus-half- 
advanced" potential). This effective potential was devised by Dirac [10] in order to handle the divergencies 
encountered in calculations of the radiation-reaction force in flat spacetime. We point out, however, that in 
general the radiative potential will not be valid in curved spacetime. Evidently, this potential yields non- 
causal results: Due to the scattering off the curvature, the Green's functions of the massless fields will have 
a support not only on the light-cone itself, but also in its interior. Thus, with the advanced potential, the 
backreaction force acting on the particle at a given moment will in general depend on the particle's future 
history. Clearly, the radiative potential will inherit this non- causal behavior, indicating that in general this 
potential is not valid in curved spacetime. Although the radiative potential yields the correct results for E 
and L z in Kerr spacetime [4] , it is not clear whether the corresponding results for Q will be correct too. 

The inapplicability of the retarded potential in curved spacetime motivates us to base our approach 
directly on the retarded potential. For concreteness and for conceptual clarity, let us refer at this stage 
of the discussion to the electromagnetic radiation reaction (later we shall generalize this discussion to a 
scalar field and implement it for explicit calculations). Consider an electrically-charged particle that moves 
in curved spacetime. In principle, the backreaction force has a simple origin: It is just the Lorentz force, 
experienced by the particle due to its interaction with its own electromagnetic field. However, the field of 
a point-like particle is divergent at the particle's location, and any attempt to calculate the backreaction 
force must first address this conceptual difficulty. In flat space, the radiative potential is often used in order 
to overcome this difficulty. Unfortunately, this method is not valid in curved spacetime. A fully-relativistic 
analysis of electromagnetic radiation-reaction in curved spacetime was carried out by DeWitt and Brchmc 
[11]. This analysis is based on the concept of a world-tube (with finite radius e ) surrounding the charged 
particle. By taking the limit e — > , DeWitt and Brchmc obtained a closed expression for the electromagnetic 
backreaction force. However, this expression is not so useful for explicit calculations: It includes a "tail" 
term, given by an integral over proper time along the particle's worldlinc. It is hard to calculate this tail- 
integral explicitly - especially because one does not know the explicit form of the Green's function that 
appears in the integrand. 

Our approach to the problem is based on the following observation: Although the particle's self field 
is divergent, the modes of the field (e.g. the spin-weighted spherical harmonics for Schwarzschild gemetry, 
and the spin-weighted spheroidal harmonics for Kerr gemetry) are everywhere finite and well-defined - even 
at the limit of a point-like source. Thus, one may simply calculate the contribution of each mode to the 
radiative evolution (by a direct application of the Lorentz-force formula to the interaction of the particle 
with the mode in question), and then sum over the modes. By this decomposition into modes we achieve 
two goals: First, we remove the singularity of the self-field. 1 Second, this is the only practical way to solve 
the field equations for a Kerr background. 



1 The contribution of each mode to the radiation-reaction is guaranteed to be finite. One still may be 
concerned about the convergence of the sum over the modes. At least for E and L z , this sum is finite. See 
the discussion at the end of the paper. 
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We still have to address the issue of renormalization. It is well known, already from flat space, that in 
general the retarded backreaction force must be renormalized before it can be used for explicit radiation- 
reaction calculations. The same phenomenon occurs for charged particles in curved spacetime. [11] (Indeed, 
it was this complication which motivated Gal'tsov [4] to use the radiative potential). Fortunately, it turns out 
that in the problem considered here no such force-renormalization is needed: The only need (and justification) 
for the force-renormalization originates from the renormalization of the mass-parameter. It therefore follows 
that, from its very nature and definition, the force-renormalization must take the form 

F« ff =F? et +Hsea a , (1) 

where F" et is the retarded backreaction force, F^a is the effective (i.e. renormalized) force, a a is the 
covariant four-acceleration, and /i se is the parameter that renormalizes the particle's mass (the particle's 
"self-energy"). This argument [and Eq. (1)] holds in both flat and curved spacetimcs. 2 In the problem 
discussed here, the radiative evolution of orbits around a Kerr black hole, no external force is considered, so 
(within the framework of the adiabatic approximation) a a = and therefore F"^ = F" et . 3 

We turn now to the concrete calculations of the radiative evolution of generic orbits in Kerr spacetime. 
For simplicity, we shall focus here on (massless, minimally-coupled) scalar-field radiation reaction. The 
generalization of this calculation scheme to electromagnetic and gravitational radiation is outlined in Rcf. 
[9]. Let C denote the constant in question - E , L z or Q . All these constants may be expressed explicitly 
as functions of coordinates and (covariant) four-velocity, i.e. C — C(x 13 ,u a ) . Thus, we may write [6] 

E = —u t , L z — u v , Q = u e + cos 9 [a (1 — u t ) + sin~ Ou^] . 

(We use Boyer-Lindquist coordinates.) The force induced by a scalar field \p(x a ) on a particle with a scalar 
charge q is given by 

F a =q(ip, a +u a u f3 ip^) = q^. a ; (2) 

this is the scalar-field analog of the electromagnetic Lorentz force [4] . [At this stage, ip(x a ) may be thought 
of as any prescribed scalar field; We shall soon apply this formalism to the backreaction problem.] It is not 
difficult to show [9] that the parameters C evolve according to 

— = C a F a =qC a ^. a , (3) 

where r is the particle's proper time and C a = fi^ 1 dC/du a . 
Next, we decompose tp into its Teukolsky modes [12] : 

^ = ^2^ , ^ = Rim^St^e^-^ . (4) 

Imuj 

From Eqs. (2,3), the contribution of each mode to the backreaction force and to dC/dr is 

F i ™=qtf™ , =qC a iP l T , (5) 

where 

The contribution of each mode l,m,u) to the change in C is thus 



AC lr 



, J C a ^ dr . (6) 



2 Note the consistency of Eq. (1) with the results obtained in Ref. [11] for electromagnetic radiation 
reaction in curved spacetime; See in particular Eq. (5.23) there. 

3 This statement can be made more precise in terms of the world- tube approach of Ref. [11] . For any 
finite radius e , the self-energy term in the backreaction force is well-defined [cf. Eq. (5.23) there] - and 
vanishes for a a = . One can then safely take the limit e — > . 
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In this expression, the integrand is to be evaluated along the particle's worldlinc x a (r). 

In order to apply Eq. (6) to the backreaction problem, we take ^p lmu> to be the retarded-field modes 
associated with the charged particle itself. To determine these modes, one decomposes the particle's charge- 
density function p{x a ) into its modes according to [12] 

47r£p = ^T^(r)^ row (%^-^ , (7) 

Imuj 

where E = r 2 + a 2 cos 2 9. The radial functions Ri m ^{r) are then determined by the Teukolsky equation [12] 

(ARi mLU ) + Vi mu) (r)Ri mLU — —Ti muJ (8) 

(with retarded boundary conditions), where A = r 2 — 2Mr + a 2 , a prime denotes d/dr, and Vi mul (r) is some 
real potential. We take the charge-density function p(x a ) to be that of a point-like particle: 

p{x a ) = qV~9~ 1 j *S(x a -x a (T))dT , (9) 

where x a (r) denotes the particle's worldline. (This approximation demands that the particle's radius d 
will satisfy d « M , which we shall assume.) Note that the thus-obtained radial functions J?; TOW (r) are 
everywhere smooth (C 1 ). 4 

Finally, the long-term evolution rate (in terms of asymptotic time t) is given by 

C = ^C lmu , C lmu = At- 1 Real(AC lmu ) . (10) 

Imu 

(The imaginary part of ACi mu will cancel out upon summation, because AC; _ m _ u = AC ; * m(i) [9].) Here, 
At = t(r2) — t(ri) , where t\ and t 2 are the lower and upper integration limits in Eq. (6). For strictly 
periodic orbits, we choose n and t 2 such that At = t 2 — n is the orbital period; For the generic orbits in 
Kerr spacetime, which are not strictly periodic but rather quasi- periodic, we take At to be sufficiently large, 
so that the quasi-periodicity will yield a well-defined limit in Eq. (10) (formally speaking, for quasi-periodic 
orbits we shall be interested in the limit At — ► oo ). 

Equations (6-10) constitute a closed scheme for the computation of the (scalar-field) radiative evolution 
of all constants of the motion in Kerr spacetime. The calculation may be further simplified, however, by 
using Eq. (9) to rewrite Eq. (6) as 



AC lri 



J EpsinfldV* (11) 



(recall V^g = Esin0 ). 5 We shall now treat the cases C = E and C — L z , which are especially simple. 
Let P denote cither E or L z , and let p denote u> or m , correspondingly. From Eq. (4) and the definitions 
of C a and C one then finds (ignoring terms which just oscillate in time) 

AP lmu = ^ [ ^ Sp sm9d 4 x a . (12) 
A 4 J 

Finally, using the decompositions (4,7), and the orthogonality of the modes, one obtains 

Pi m . = g J R lm " (rW* mu , (r)dr + C.C. . (13) 



4 The special case of circular orbits (r=const) is exceptional, as the radial functions are not C 1 at the 
particle's location. This simpler case has to be treated separately. 

5 C a is only defined along the particle's worldline; However, the integrand in Eq. (11) vanishes elsewhere 
anyway. 
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It would be interesting to compare Eq. (13) to the standard expressions for energy and angular- 
momentum radiative losses. By virtue of Eq. (8), it is possible to integrate Eq. (13) explicitly: 

Pi mu = g [Afl* 1 ™^;™]™ +C.C. , (14) 

where r + denotes the horizon's radius (this equation may easily be checked by a direct differentiation). 
Taking into account the asymptotic behavior of the radial functions in both asymptotic limits [7] , it is now 
straightforward to show that Eq. (14) is exactly of the desired form; Namely, Pimui is the sum of two terms 
- the fluxes (per unit mass p ) of energy or azimuthal angular momentum that flow (i) to infinity and (ii) 
into the horizon - exactly as calculated by the standard method [7] . 

The more interesting application of our scheme would be the calculation of Q (here there are no previous 
results to compare with). Equations (10-11) provide a closed scheme for the calculation of Q . In Ref. [9], 
however, we derive from Eqs. (10-11) a more explicit expression for Qi muJ : 

Q lmuJ = ^ J (HR lm » Tf mu - R lm " ft mio - R l ™ft mio ) dr + C.C. , (15) 

where 

H = irA^ 1 [ujE{r 3 + ra 2 + 2Ma 2 ) - 2(loL z + mE)Ma + mL z (2M - r)] 

and Timu (r) and Ti muJ (r) are functions analogous to Ti mu , obtained by replacing the factor p at the left-hand 
side of Eq. (7) with p = 2ru r p or p = Au r p , correspondingly. Equation (15) thus allows a straightforward 
numerical calculation of Qi mu> (and Q ). 

We wish to emphasize that no attempt was made here to provide a complete and fully-rigorous treatment 
of all the subtle issues involved in the radiation-reaction problem. Rather, the goal here was to propose a 
simple and straightforward approach to the problem, which appears to yield promising results: It recovers 
the standard results for E and L z in Kerr spacetime (at least for a scalar field) - and also yields an explicit 
expression for Q . Our approach is based on the delta- function approximation, Eq. (9). This approximation 
has extensively been used over the years for radiation-reaction calculations (see for example Refs.[3-5] and 
references therein). We believe that this approximation is justified in our case. We hope that in the future 
it will be possible to apply more rigorous methods for the calculation of the radiative evolution of orbits in 
Kerr spacetime. 

One may be concerned about the convergence of the sum over the modes in Eq. [10] . This series is 
known to converge for C = E and C — L z , and this makes us optimistic about its convergence for C = Q 
as well. The explicit calculation scheme presented here allows one to verify this convergence (for Q = C) 
numerically - and perhaps even analytically. 

Finally, we mention several tests which may be applied to our result for Q , Eq. (15). The first such 
test would be, of course, to verify converges of the sum over the modes (see above). Another simple test can 
be applied to the special case a=0 (the Schwarzschild case): Due to the spherical symmetry of the latter, 
the orbital plane must be preserved. A simple calculation then yields Q = 2(Q/L Z ) L z . Although this test 
is limited to the case a=0 , it is conceivable that if our scheme is faulty, it will yield wrong results in that 
case too. Finally, for equatorial orbits (Q = 0) in Kerr spacetime, Q must vanish identically 

I would like to thank Eric Poisson, Sam Finn, Eanna Flanagan and Kip Thorne for stimulating discus- 
sions and helpful comments. 
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